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Signatures of entanglement in an optical tomogram 
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We theoretically study the optical tomography of maximally entangled states generated at 
the output modes of a beam splitter. We consider even and odd coherent states in one of the 
input modes and vacuum state in the other input mode of the beam splitter. We have shown 
that the signatures of entanglement can be observed directly in the optical tomogram of the 
state, without reconstructing the density matrix of the system. Two distinct types of optical 
tomograms are observed in any one of the output modes of the beam splitter based on the 
quadrature measurement in the other output mode if the output modes are entangled. The 
different features shown by the optical tomograms are verified by investigating the photon 
number statistics of the corresponding state. 
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I. INTRODUCTION 


Quantum entanglement is a key resource for quantum information and quantum computing. 
After the EPR paper 111, a tremendous amount of work has been done in the field of quantum 

n 

entanglement [2|. In most of the quantum information processes, such as, quantum teleportation 
jjj, quantum cryptography jj], quantum metrology jfj], etc., the systems are prepared initially in 
an entangled state. Much attention is devoted to the discussion of entanglement properties of 
continuous variable systems, for their great practical relevance in applications to quantum optics 
and quantum information. Generating entangled states is one of the major tasks in quantum 
information processing. Various devices have been proposed and realised experimentally to generate 
quantum entanglement. The simplest one is by using a quantum mechanical beam splitter. 

Once the entangled states are created in an experiment, it is important to characterise the 
state of the system precisely. Optical homodyne tomography can serve as an efficient technique 
to measure and reconstruct the state of entangled optical fields. Optical tomography is based on 
the one-to-one correspondence between the quasiprobability distribution and probability distribu¬ 
tion of the rotated quadrature phases of the optical field [6j. The first experimental observation 
of squeezed state of light, by measuring the quadrature amplitude distribution using balanced 
homodyne detection, has been done in 7|. A series of homodyne measurements of the rotated 
quadrature operator of the field develops a quadrature histogram called optical tomogram. The 
optical tomogram contains complete information about the system and can be considered as a pri¬ 


mary object which describes t 
density matrix of the system 


re system completely, other than the conventional state vector or the 
]. After this experiment [j|, many nonclassical states of light have 
been characterised using optical homodyne tomography. A review of continuous-variable optical 
quantum-state tomography is given in 9|. 

Entangled states of light have been recently[lCHl3l] characterised using optical homodyne tomog¬ 


raphy. In these experiments, the two-mode density matrix of the system has been reconstructed 
from the optical tomogram and the entanglement is quantified using the reconstructed density 
matrix. A conditional measurement on one of the modes of entangled states may change the state 
in the other mode due to entanglement, and such changes may shows up in the optical tomogram 
of the state. In this paper, we study the optical tomography of bipartite continuous-variable en¬ 
tangled states of light. Main goal of this work is to find the signatures of entanglement in the 
optical tomogram of the state, without reconstructing the density matrix of the state. For this 
purpose, we investigate the optical tomogram of maximally entangled coherent states created in a 
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beam splitter. The analytical results presented in this paper will be useful for the characterisation 
of the continuous-variable entangled states. 

The optical tomogram of a quantum state can be theoretically evaluated using the symplectic 
tomograhic approach 14- 19f]. The optical tomogram calculated by this method can be used to 


compare and verify the optical tomogram obtained by the homodyne detection of photonic states. 
The symplectic tomography of several single-mode nonclassical states of light have been investi¬ 


gated in the literature 


10-23]. A theoretical investigation of the optical tomogram for two-mode 
coherent states of charged particle moving in a varying magnetic field is available in Q . The cor¬ 
rectness of measured tomogram can be checked using the properties of tomogram, like uncertainty 
relations [251 ]. tomographic entropic inequalities [261 ]. purity constraints [27], etc. An experimental 


check of the two-nrode Robertson uncertainty relations and inequalities for the highest quadrature 
moments using homodyne photon detection have also been suggested [28|. The operational use of 


experimentally measured optical tomograms to determine state characteristics avoiding the recon¬ 
struction of quasiprobabilities has been demonstrated in [ 29 ]. We use the symplectic tomographic 
approach to find the optical tomogram of the two-nrode photonic state. This paper is organized 
as follows. In section hd we discuss the generation of maximally entangled states using a beam 
splitter. Section IIIII is devoted to the discussion of optical tomography of entangled states gen¬ 
erated at the output of the beam splitter. It also explains the signatures of the entanglement in 
the optical tomogram of the state. Section IIVI gives a verification for the different features shown 
by the optical tomogram by investigating the photon number statistics of the corresponding state. 
Finally in section [3 we summarize the main results of this paper. 


II. BEAM SPLITTING ACTION 

Consider a 50/50 beam splitter with zero phase difference between reflected and transmitted 
port. The unitary operator for the beam splitter reads 


U = exp 


— {a)b — ab t) 

14 


(1) 


where a and b are the bosonic operators for the input field modes. The output field modes of 
the beam splitter are designated by c and d. A schematic diagram of the beam splitter is given 
in fig. [U A beam splitter generates entangled state if one of the input fields is nonclassical H- 
We consider both classical and nonclassical states in the horizontal input port (mode a) of the 
beam splitter and study the optical tomogram of the output states. In both of these cases we take 
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FIG. 1. A 50/50 beam splitter with \ip) in the horizontal port and |0) in the vertical port, a and b (c and 
d ) are the input (output) field modes. 


vacuum state |0) in the vertical input port (mode b ) of the beam splitter. The states considered for 
classical and nonclassical states are coherent state, and even and odd coherent states, respectively. 
The initial even and odd coherent states exhibits different nonclassical behaviour during the time 
evolution when compared to an initial coherent state [32J]. Beam splitting action on the coherent 
state |V 7 ) = |a), where a (= |a| e lS ) is a complex number, with vacuum |0) generates the separable 
state 


l*>„ = \P) ® \P) , (2) 

where /3 = a/\/2. Next, we consider beam splitting of even and odd coherent states, defined by 


= N h 


\a)+e inh \-a) 


where the normalisation constant 


1 




l + (-l) h e" 2|Q|2 


( 3 ) 


( 4 ) 


Here h = 0 and 1 corresponds to the even and odd coherent states, respectively. In this case, we 
get entangled states at the output modes of the beam splitter. The state |4>) /) of the beam splitter 
output modes is calculated using the unitary operator given in eq. ©: 


m h = N » 


+ e inh \—0) c \ —@) d 


( 5 ) 


In the limit of large |a| , the coherent states appearing in the superposition given in eq. @ form 
an orthogonal basis and thus the entangled state given in eq. (|SJ) can be written in the Schmidt 


decomposition form 


31]. Since all the Schmidt coefficients have the same magnitude l/\/2, the 
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entanglement of the state 1$}^ can be found to be E = log 2 (2) = 1 ebits, which is the maximum 
entanglement possible in 2 dimension. Hence, for large |a| 2 , the state is a maximally entangled 
state in 2 dimension. 


III. OPTICAL TOMOGRAPHY 


We theoretically calculate the optical tomogram of the maximally entangled state \&) h and sep¬ 
arable state |d>) cs using the symplectic tomogram of the state. Consider the homodyne quadratures 
X and X^ 2 ,v 2 associated with mode c and mode d respectively: 


X^j^ = iii qi + zq f)\, X ll2[U2 = p 2 q 2 + ^2 P 2 - 


( 6 ) 


Here pi, zq, p 2 , zq are real parameters and qs and ps are the usual position and momentum 
operators. The symplectic tomogram of a two-mode state with density matrix p is defined as 

m 


Af (X ^, pi , zq, X ^ 2 ^ 2 , p2 j ^2) 
= Tr 


P&( x lii,v\I X^ 1Ul 


) ${Xfj, 2t v 2 1 Xfj, 2t i/ 2/ 


( 7 ) 


where X )M , A and X^ 2yU2 are the eigenvalue of the Hermitian operators X fn Ul and X^ 2 ^ V2 , respec¬ 
tively. For a pure two-mode state with wave vector |<I>}, the symplectic tomogram can be obtained 
by taking the fractional Fourier transform of its coordinate wave function <3?(x,z/) [24, 33|: 


Xd (A /^1 ^1 , , Zq , X^ 2 p 2 , P2 ; ^2 ) j __ 2 


1 


/ 


TW 2 I 

dx dy <b(x, y ) exp ( ^-x 2 — ^2£biXL x 
2 zq zq 


( 8 ) 


ip 2 2 iXf i2 v 2 

x ex P 7^— v --—y 

2^2 


In experiments optical tomogram of the states are generated and we look for the signatures of 
entanglement in the optical tomogram. The optical tomogram of a two-nrode state can be obtained 


from the symplectic tomogram given in eq. j8j) by the following substitutions 


nay, 


33]: 


-k //1 A/y — Xq 1 , pi — COS01, v 1 — sin0i, 

X H 2 ,V 2 = x e 2 , P 2 = cos 02, zq = sin 0 2 . (9) 

Here Xq 1 and 9\ ( Xg 2 and 9 2 ) are the quadrature and the phase of local oscillator in homodyne 
detection setup for mode c (mode d), respectively. The phase of the local oscillators varies in the 
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domain 0 < 9\, 62 < 2n. The optical tomogram of the two-mode state is then reads as 


u(X ei ,0 1 -X 62 ,d 2 ) 

= M ( Xq 1 , cos 9 1 , sin 9\ ; Xq 2 , cos 62 , sin @ 2 ) • 


( 10 ) 


Two-mode coordinate wave function <f>h(x,y) for the entangled state |4>) h is calculated as 


* h (x,y) = ?±J2e- i * rh x 
v 71 n 

r=0 


exp(— |/3| 2 /2 — /3 2 /2 + V2/3 r x — x 2 /2)x 
exp(— |/3| 2 /2 - /3 2 /2 + \/2/? r y - y 2 / 2), 


( 11 ) 


where /3 r = /3e l7Tr . Substituting &h(x,y) in eq. (JSJ) and using the transformation given in eq. (|1U|) . 
the two-mode optical tomogram for the entangled state is obtained as 


where 


^h (-^Q?i; ^11 X62 > ^2) 


7 r 


£■ 

r=0 


—inrh 


Vr(Xg 1 , 9i)r}r {Xg 2 , 62 ) 


rir(X 0 i ,9i ) = exp 




+\Z2/3 r X 0i e-^ 



( 12 ) 


( 13 ) 


for i = 1 and 2. 

Now, for the separable two-mode state |<£) cs = |/3L |/3} d the optical tomogram can be written as 
the product of optical tomograms of the subsystems [8l], that is u cs (Xg 1 ,9 \; Xg 2 ,02) = uq (X^, 9\ ) x 
uj 2 {X() 2 . 92 ); where uii(Xg 1 ,0i), uj 2 {Xg 2 : 92 ) are the optical tomogram for mode c and mode d, 
respectively. It can be shown that 


u i {Xg i ,9 i ) = —= exp 

' 7T 


- (x fl .-x/2|/3|cos(<J-0i))‘ 


( 14 ) 


where i = 1 or 2. The maximum intensity of the optical tomogram oj^Xg.. 9 t ) is l/y/n, which 
occurs along a sinusoidal path, defined by Xg i = \J 2 |/3| 2 cos(5 — 9i), in the Xg. — 9i plane. Hence 
the optical tomogram in mode c (mode d ) will always be a sinusoidal single-stranded structure for 
any values of parameters in mode d (mode c) and 6 . 
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A. Signature of entanglement 


The state \&) h is an entangled state and the corresponding optical tomogram cannot be written 
in the form of the product of subsystem tomograms, that is uJh{Xg 1 .0 \; Xq 2 . 02 ) / (jJ\{Xq 1 ,6\) x 
u) 2 (Xg 2 , O 2 ). We present our analysis for h = 0, that is, for initial even coherent state. When h = 0, 
the entangled output state |<h) 0 = No [|/3) |/3) + |— /3) | —/3)]. A measurement of field quadrature in 
any of the modes will collapse the entanglement between the modes. That is, the quantum state 
in one mode is correlated with the quadrature measurement in the other mode. For example, a 
measurement of quadrature Xg 2 in mode d will project the state |d>) 0 to the state 1 4>) 0c in mode c: 

l<A)o,c = ^0 [M x e 2 M |/3) +ll>-f,(Xg a , 6 2 ) |-/3)], (15) 


where No is appropriate normalisation constant and ip±p(Xe 2 , 6 2 ) = (Xg 2 ,02 |d=/3), is the quadra¬ 
ture representation of the coherent state |±/3). Based on the relative strength of the coefficients 
ipp(Xg 2 ,d 2 ) and ip-p(XQ 2 , 62 ), the state |</>) 0 c can be one of following: |/3), \—/3) and a superposition 
of |/3) and | —/3). For large |a| 2 , the probability for occurring the state |±/3) is proportional to 


\^±p(Xg 2 ,02 )\ 2 = —j= exp 
V7T 


|/3| 2 — 2 |/3| 2 cos 2(6 — 62 ) 


-X $ 2 ± 2V2X d2 |/?| cos (5 - 0 2 ) 


(16) 


For Xg 2 / 0, relative strength of the probabilities crucially depend on the last term in eq. (fl6l) . 
In the range n/2 < \5 — 62 ] < 37r/2, the state \4>) 0c becomes a coherent state |—/3) because 
\'i/j—p(Xg 2 ,^ 2 )| 2 \' 4 >/ 3 (Xq 2 , 02 )\ 2 ■ Thus, the optical tomogram in mode c will be a single-stranded 

structure corresponds to the coherent state \—/3) (Note that the optical tomogram for a coherent 
state |/3) in Xq. — 0 t plane is a sinusoidal single-stranded structure as given in eq. (fTTl) ). Also, 
in the range 0 < \5 — 62 ] < vr/2 and 37r/2 < |5 — 62 ] < 2vr, the state | (fi) 0c becomes a coherent 
state |/3) because \tpp(Xg 2 , $ 2 )| 2 \ip-p{Xg 2 , @ 2 )\ 2 , which gives a single-stranded structure for the 
optical tomogram in mode c. Figure [2]( a) shows single-stranded structure in the optical tomogram 
ojq(Xq 1 ,0\\ Xg 2 , 62 ) for |a| 2 = 10, A = 0.2, Xq 2 = 2.0 and |<5 — 021 = 0.3. 

It can be shown that \^^{Xq 2 , ^ 2 ) | 2 = \fp~p{Xg 2 , ^ 2 )| 2 for |<5 — 021 = 7 ^/2 or 37r/2 (within the 
periodicity of optical tomogram). For these two values of 1 5 — 021, the probability for occurring |/3) 
and 1-/5) in mode c is 50 : 50 and hence the state |</>) reduces to even coherent state of the form 
[|/3) + |—/3)]. The optical tomogram in mode c will display double-stranded structure, in which, 
one strand corresponds to |/3) and the other corresponds to | —/3). The optical tomogram in mode c 
for |<5 — 021 = 7 t/ 2 and 37 t/ 2 with Xg 2 = 2.0 are shown in figs. EXb) and[2](c), respectively. Quantum 
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FIG. 2. Optical tomograms in mode c for the entangled state |$) 0 with \a\ 2 = 10, 

6 = 0.2, and Xg 2 = 2.0, for different relative phases |<5 — 0 2 \ of the quadrature measurement in mode 
d: (a) 0.3 (b) n/2 and (c) 37r/2, respectively. Optical tomogram is a sinusoidal single-stranded structure 
for \5 — 0 2 \ = 0.3. For \S — 0 2 \ = 7r/2 and 37t/2, the optical tomogram shows sinusoidal double-stranded 
structure. 

interference between the state | j3) and \—/3) are reflected in the optical tomogram at regions in 
Xg 1 — 61 plane, where the two strands intersect. When Xq 2 = 0, the state |<^) 0 c is always an even 
coherent state, without any condition on 15 — ^ 2 1 - This displays double-stranded structure in the 
optical tomogram. 


IV. MANDEL’S Q PARAMETER 


We can also show the above features using the statistics of 


loton number in the state 


specifically, in terms of the Mandel’s Q parameter, defined as 

q = (fl)2 - 1 , 

(n) 


M 


'0,c> 


(17) 


where h is the photon number operator. A positive value of Q indicates the super-Poissonian 
statistics of the field, and Q = 0 indicates the Poissonian statistics exhibited by a coherent field. Q 
parameter of the state |</>) 0 c in mode c as a function of relative phase \5 — 62 ], is plotted in fig. [3j 
It shows that at |<5 — @ 2 ! = 7 r/ 2 , 3ir/2 and close to it, the state exhibits super-Poissonian statistics 
and for all other |<5 — 6 * 2 1 values, the state |</>} 0c shows Poissonian statistics corresponding to a 
coherent field, which is either |/3) or |— j3). The different statistics of photon number exhibited by 
the state in mode c upon changing the parameters in mode d, has been experimentally observed in 
the case of micro-macro entanglement of light jp]] using the reconstructed density matrix. We have 
theoretically shown that, without reconstructing the density matrix of the system, the signature 
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|6 - 6 2 1 


FIG. 3. Mandel’s Q parameter for the state 1 4>) 0c as a function of the relative phase |<5 — ^ 2 1 of the 
measurement in mode d with \a \ 2 = 10. The positive value of Q at |<5 — d 2 \ = 7t/2 and 3n/2 (and close to 
it) indicates the super-Poissonian statistics of the state |^) 0 and for all other values of |<5 — $ 2 !, the state 
|</>) 0 c exhibit Poissonian statistics (Q = 0). 



FIG. 4. Optical tomograms u>i(Xg 1 , $i; Xg 2 , 0 2 ) hi mode c for the entangled state |<f>) 1 with \a \ 2 = 10, 
6 = 0.2, and Xg 2 = 2.0, for different relative phases |<5 — 021 of the quadrature measurement in mode 
d: (a) 0.3 (b) n/2 and (c) 37r/2, respectively. Optical tomogram is a sinusoidal single-stranded structure 
for \5 — 0 2 \ = 0.3. For \S — d 2 \ = n/2 and 3n/2, the optical tomogram shows sinusoidal double-stranded 
structure. 


of entanglement can be directly observed in the optical tomogram of the state. When the initial 
state is an odd coherent state (i.e., h = 1), we get the entangles state | < h) 1 given in eq. (HJ at 
the output modes of the beam splitter. We have repeated the forgoing analysis for the entangled 
state |<h) 1 and verified the results obtained earlier. The figs.[I](a)-(c) shows the optical tomograms 
ui(Xg 1 , 6 \; Xg 2 , 62 ) in mode c for the entangled state | < h) 1 with same set of parameters used in the 
case of the entangled state |d>) 0 . 
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V. CONCLUSION 

A closed form analytical expression for the optical tomogram of the maximally entangled coher¬ 
ent state generated at the output of the beam splitter is derived by taking the fractional Fourier 
transform of the output two-mode wave function. For separable two-mode states, the optical to¬ 
mogram of the system can be written as the product of the optical tomograms of the subsystems. 
Whereas, for the entangled two-mode states, the optical tomogram in the mode c shows different 
features when we change the parameters Xq 2 and 62 in the mode d. Similarly, the optical tomo¬ 
gram in the mode d will be affected by the parameters in mode c. Specifically, for the entangled 
state |<h) h = Nh [|/3) |/3) + e mh |—/3) |—/3}], with Xg 2 / 0, the optical tomogram in mode c shows 
double-stranded structure if |<5 — 6*2 1 = 7t/2 or 3ir/2 and a single-stranded structure for all other 
values of |<5 — 62 |. Our analytical result for the optical tomogram of entangled coherent states may 
be useful for the experimental characterization of such states. 
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